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Abstract 

We investigate behavior of the Fisher information matrix of general stable distributions. Du- 
Mouchel (1975, 1983) proved that the Fisher information I aa of characteristic exponent a diverges 
to infinity as a approaches 2. Nagaev and Shkol'nik (1988) made more detailed analysis of I aa and 
derived asymptotic behavior of I aa diverging to infinity as a approaches 2 in the symmetric case. 
Extending their work in this paper we have obtained behavior of the Fisher information matrix of 
general stable distributions as a approaches 2 by detailed study of behavior of the corresponding 
density and its score functions. We clarify the limiting values of the 4x4 Fisher Information 
matrix with respect to the location //, the scale cr, the characteristic exponent a and the skewncss 
parameter (3. 
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1 Notations and preliminary results 

The family of stable distributions has enjoyed great interest of researchers in many fields e.g., math- 
ematics, physics, cosmology and even economics. These applications are summarized in Uchaikin 
and Zolotarev (1999). In statistical inference estimation of stable parameters has been of great inter- 
est. In recent years, maximum likelihood estimation of stable distributions has become feasible (see 
Brorsen ans Yang (1990), Nolan (2001) or Matsui and Takemura (2004)). Even in time series models 
like GARCH using general stable distributions, maximum likelihood estimation is possible owing to 
recent development of global algorithm (see Liu and Brorsen (1995)). Since the Fisher information 
matrix gives useful criteria for the accuracy of estimation, it is indispensable to analyze that of general 
stable distributions. Near Gaussian distribution (a = 2), the information of a, I aa diverges to oo 
and asymptotic behavior of I aa as a f 2 is of great interest. Nagaev and Shkol'nik (1988) have solved 
this problem excellently for symmetric stable distributions. However for general stable distributions 
Nagaev and Shkol'nik (1988) stated "We note the problems under study are as yet unresolved for 
non-symmetric stable distributions." In this paper we investigate and obtain asymptotic behavior of 
lee, = /j,,a,a, (3, as a f 2 under general stable distributions. 
Let 

$(i) = $(t;(j,,a,a,(3) = exp (-\at\ a j 1 + i/3 (sgn t) tan (J^-j (|o^| 1_a - 1)} + ifd 
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denote the characteristic function of general stable distribution (a/1) with parameters 

= $!, 62,63,64) = (n,a, a, 0), 

[i e R, a > 0, < a < 2, \0\<1. 
For the standard case (fi,a) = (0,1) we simply write the characteristic function as 

(1.1) $(t; a, /3) = exp {l + i(3(sgnt) tan (^) (I*! 1 -" - 1)}) . 

This is Zolotarev's (M) parameterization (see p. 11 of Zolotarev (1986)). The corresponding density 
is written as f(x; fi, a, a, /3) and f(x; a, f3) in the standard case. Then 

f(x;n,cr,a,f3) = -f ( - - ;a,PJ . 



1 / x 2 



We also write the density of N(Q, 2) as 

/(»;2) = 77^ exp 

For the rest of this paper, without loss of generality we consider the Fisher information only at the 
standard case (//, a) = (0, 1). 

The first derivative of f(x; a, ft) with respect to x is denoted by f'(x; a, f3) and the partial deriva- 
tive with respect to &i, i.e., 

df(x;n,a,a,P) 



(/i,a)=(0,l) 



d6i 

is denoted by fo^x; a, (3). Note that 

(1-2) f„(x;a,(3) = -f'(x;a,(3), 



(1.3) f a (x; a, (3) = -f(x; a, f3) - xf'(x; a, (3). 

The components of Fisher information matrix / are defined as follows. 

[°° df df 1 



d9 l 86 j f 



Let 

/ir a 



C = —(3 tan 



2 / fna 

n = — arctan p tan — 
H ira Y V 2 

From Theorem 1 of Nolan (1997), which is a modified version of the formula (2.2.18) of Zolotarev 
(1986), the density f(x; a, /?) for the case a ^ 1 and x 7^ £ is written as 

| x _£|l/(a-l) /-i 

-0 



(1.4) ^ g ' a '^ = 2la-ll / ^a.^exp^lx-Cr^-^A^a.jS))^, 
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where g* = gsgn(x — £) and 

cos(f^) \~ cos + (a-l)<p}] 



i 

7T \ \ T^T 



(1.5) A((p;a,P) = (cos I -ag , , ■ . f7r 

is a positive function. Here we should note the relation f(x;a,f3) = f(—x;a,—f3) for x — ( < 0, 
which will be used many times. As a f 2, £ and g converge to and at a = 2 we obtain an unusual 
representation of Gaussian distribution N(0, 2), 

f(x;2)=x [ 1/ (2 sin exp (-x 2 / (2 sin ^ x > 0. 



The remainder of this paper consists of three sections. In Section 2 we derive behavior of 
f'(x;a,{3), fp(x; a, (3) and f a (x;a) as a | 2 by careful analysis of the formula p.4|) . In Section 
3 behavior of f a (x;a,(3) and fp(x;a,0) are obtained by utilizing the inversion formula. The infor- 
mation matrix of general stable distributions is given in Section 4. 

2 Density and derivatives of general stable distributions close to 
the normal distribution 

In this Section the density and derivatives of general stable distributions as a f 2 and x — > 00 are 
obtained. In the following we write 

A = 2- a. 

As A — > we can show 

(2.1) f(x;a,(3) = f(x;2)+0(A) 

uniformly in x by the finiteness of f a (x; a, (3), which can be easily verified by the formula (|3.2|) in the 
proof of Lemma 13.11 However if we consider the ratio /(x; 2)/f(x; a, (3) rather than the difference 
fix; 2) — f(x; a, (3) of the densities as x — > 00, f(x, 2) is much smaller than f(x; a, /?) and we can not 
obtain an accurate approximation of f(x,a,f3) without closer investigation of the remainder 0(A). 
In the case of x — > 00 we can utilize the asymptotic expansion 

(2.2) f(x; a, 0) = A(l + (3){x - C) A " 3 + o(Ax" 4 ) 

from the expansion (2.5.4) on p. 94 of Zolotarev (1986), Bergstrdm (1953) or Section XVII. 6 of Feller 
(1971). (2.5.4) of Zolotarev (1986) is in terms of (B) representation with characteristic function 
defined on p. 12 of Zolotarev (1986). Hence as x — > 00 

f B (x; a, 13) = A(l + (3 B )x A ~ 3 + o(Ax~ 4 ), 

where fs and (3b means the density or (3 of (B) representation. We obtain (|2.2|) by using the relation 

f(x;a,f3) = (cos (^ap)) " !b ^(cos (t^)) " ( x ~ ();u,Pb 

by noting that cos^a^) — ► 1 and 0b = P + 0(A 2 ) as A — > 0. Then the problem is which of the two 
approximations (|2.1|) and (|2.2j) is dominant for large x. We give an answer to this problem in the 
following two theorems, which are based on Theorem 1 of Nagaev and Shkol'nik (1988). We can also 
see the precise statement on p. 129 of Uchaikin and Zolotarev (2003) concerning this problem. 
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Theorem 2.1 Let (3, — 1 < f3 < 1, be fixed and let (3* = (3sgn{x — £). Define 

F 1 (x;a,P) = f(x-(;2), 
F 2 (x;a,P) = (l+/T)A(x-C) A ~ 3 , 
g(x;a,(3) = Fi(x; a, (3) + F 2 (x; a, (3). 

For an arbitrarily small e > there exist Aq and xq such that for all A < Ao and \x\ > xq, 

\f(x;a,P)/g(x;a,/3) - 1| < e. 

Furthermore, for an arbitrarily small constant 5 > 0, 

F 1 (x;a,f3)(l + o(A 5 / 2 )) if \x - (\ < (2 - <5)(log 1/A) 1 / 2 , 

F 2 (x; a, (3) (l + o(A 5 / 2 )) if |x - C| > (2 + <5)(log 1/A) 1 / 2 , 



g(x;a,/3) 



uniformly in \x\ > xq. 
Theorem 2.2 Define 

F[(x;a,(3) = f(x - C; 2), 

^(x;a,/3) = _3(l + /?*)A(x-C) A - 4 , 
g'(x;a,f3) = F[(x;a, (3) + F 2 {x;a, f3). 

Under the same conditions and notations of Theorem \2.1l for an arbitrarily small e > there exist 
Aq and xq such that for all A < Ao and \x\ > xq, 

\f'(x;a,(3)/g'(x;a,(3) - l| < e. 

Furthermore, for an arbitrarily small constant 5 > 0, 

' F[(x;a,P) (1 + o(A< 5 / 2 )) if \x - C| < (2 - <5)(log 1/A) 1 / 2 , 

< 

k ^(x;a,/3)(l + o(A< 5 / 2 )) if \x - C| > (2 + 5) (log 1/A) 1 / 2 , 



g'(x;a,(3) 



uniformly in \x\ > 

We need the following 4 lemmas concerning behavior of ^4(</9; a, /3) in (jl.5|) for ip = 1 to prove 
Theorem 12 . 1 1 and Theorem l2.21 which correspond to the lemmas of Theorem 1 of Nagaev and Shkol'nik 
(1988). In the lemmas and the proofs of Theorems 12 . 1 1 and 12 . 2 1 e and e' denote arbitrarily small positive 
constants because in the end of the proofs of Theorems 12 . II and 12 . 21 we let e and e' j separately. From 
now on the Fisher information is evaluated at a fixed value of f3, \f3\ < 1. For notational abbreviation 
we sometimes write A(ip;a, (3) as A(ip). We also write 

A = 1 - <p, <pa = 1- A 1/2 ~ e . 
Lemma 2.1 As A -> and for < A < A/e' , 

uniformly in A. 
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Lemma 2.2 isA^O 

A(<pa) = \ + ^A 1 ^ + o(A), A'( VA ) = -yA 1/2 " £ " ^±^AV^e + 0(A) _ 
Lemma 2.3 isA^O and A ->■ swc/i i/iai A > A 1 / 2_e , 

v , 1 7T 2 2 n2 . 7T 2 , A 2 (l-A + /?) 2 /x2 n 

^4(1 - A) = - + —A 2 + o(A 2 ), A'(l — A) = — —A v gA3 ^ + o(A 2 ), 

4"n n- f2 / 2 i2 3 A 2 (1-A + /?) 2 3 A 2 (1-A + /?) / A(l-A + /3) 2 ^ 2 
A (l-A)-- + y A Xi - 3 ^1 - 2 j+o(A ). 

Note that A/e' < A 1//2_e < A implies A = o(A 2 ) and we have to consider terms like A 2 /A 4 or A 2 /A 3 
in Lemma EH 



Lemma 2.4 If A is sufficiently small, A(ip) is a monotonically decreasing function on ( — g*,l). 

The proofs of the lemmas are given in Appendix. 

In the following c > is an appropriate positive constant. 

Proof of Theorem 2.1 We prove the assertion in the case x — C > first. As in the proof of 
Theorem 1 of Nagaev and Shkol'nik (1988), let t,t] > be arbitrarily small numbers. We denote 
z = (x — £) a /( a_1 ) and tpo = ipA — z~ l / 2+T . Then we divide the integral 

H= / A(ip;a, (3) ex.p(—zA(ip;a, (3))dip 

J-Q* 

into six subintegrals, 

6 
k=l 

where each Hf~ corresponds to the integration of H for the k-th. interval of [-Q*, 1 — rf), [1 — r\, ip®), 
Ypoi Va)) [VAi 1 — A/e'), [1 — A/e', 1 — Ae') and [1 — Ae', 1]. First we calculate H[s utilizing Lemmas 
l2.1ll2.4l for fixed e and e' . Later we see that H% and Hq dominate the others terms. We also let e and 
e' I in the end. 

Calculation of Hi: 

From Lemma 12.21 A(ipa) > \- By Lemma 12.41 and 1 — rj < i/ja, for sufficiently small A there exists a 
constant pi £ (0, 1) such that piA(l — if) > j. Then for a constant 7 > \ we easily find 

(2.3) Hi = — ^ rf \(l-pi)A{ip)exv{-zpiA(<p)-z{l- P i)A{ip)}dip 

z{l - pi) J- g * 

< — r exp{-zpiA(l - rf)) 

2(1 - Pi) 
= 0(exp(— r yz)/z). 
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Calculation of H^- 

From Lemma 12. 



Then we can write 



r<po 

H 2 = A(ip)exv(-zA(tp))dip 

Jl-T) 

< r]A(l - T])exp(-zA((f )) 
= 0(exp(-zA(<po))). 



A{w) = A(<p A ) + A'(^ A ){ipo - <p A ) + -A"(£)(<po - v? A ) 2 , 



where tpo < £ < <p A . By Lemma 12.21 and Lemma 12.31 we can write 

zA{<p A ) = | + i?i(A,z,e), 

zA'((p A )((po ~ Va) = #2(A,2,e), 
zA"(0(W)-VA) 2 < cz 2t , 

where 

fli(A,z,e) = O(zA 1_2e )>0, 
2fc(A,s,e) = 0(A 1 / 2_e z 1 / 2+T ) > 0. 

Here we use functions Ri(A, z,e), i = 1, 2, for convenience in comparing #2 with #3. Then 
(2.4) J ff 2 = o(exp(-^- J Ri(A,z,e)- J R 2 (A,z,e)-cz 2r )) . 

Calculation of H%: 

We show the detailed calculation of #3 because this is the dominant term in H as A — > and x — > 00. 
We have for £ £ [i^o, <£>a], 

#3 = exp(-zA(ip A )) / ^(v 3 ) ex P ( -zA'(ip A )((p - ip A ) — -A"(^)((p - (p A ) 2 ) dip. 
Utilizing Lemmas 12.21 and 12.31 for some constants p\,p2 £ (0, 1) we obtain 

2 

-zA(ip A ) = -- - ^-A 1 ~ 2e z + o(Az), 
4 lb 

-zA'(cp A )(ip - ip A ) = - 



-A 1 / 2 - + -(1 + /3) 2 A 1 /2+3^ piZ V2+T + (Az 1 / 2+T ), 
8 8 / 

A"(0(<P-<PA) 2 = -^y(^-^o) 2 

/vr 2 2 3 A 2 (1-A + /?) 2 3 A 2 (l-/? + A) / Ajb-Uf \ 2 \ p 2 2t 
"^Y A "8 A^ 4 A^ V A*" ) +o{X) JY Z ■ 
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Summarizing the main terms, we obtain 



A(<p) = \ + 0(X 2 ), 
■zA(<pa) = ~-Ri(A,z,e) 



-zA'((p A )((p - (p A ) = -R 3 (A,z,e), 

- Z -A"(^-^a) 2 = -|^(^-^o) 2 + 0(A 2 ^) + 0(A 2 /A 4 z 



^(^-^) 2 + 0(z-^) + 0(AV), 



where 



< R 3 {A,z,e) = 0(A 1 / 2 ~ e z 1 / 2+T ) < R 2 (A,z,e). 
Note tp < ip < ip a & A 1 / 2 ' 6 < A < A 1 / 2 "' + z~ 1 / 2+T '. Then we obtain 



H 3 = R(X, A,z,e)^exp {-^j j exp y 



Z 7T^ 

■^(v-^a) 2 ) dtp, 



where 

R(X, A, z, e) = (1 + 0(A 2 )) exp (-i?i(A, z, e) - R 3 (A, z, e) + 0(z" 1+4r ) + 0{A 4t z 2r )) . 
Here the inequality 

POO 1 

/ e ~y 2 n dy < i e ^ 2 / 2 , x>0 

gives 

Vo eX H"2T ( ^-^ o) ) dip = TJo eXP l"2^ ) dip -TzL GX H"2^ 

2 _/ 1 (-z 2T 

+ °( 7v^ exp 



7TC 



Then we obtain 

(2.5) F 3 = exp (-| - ^(A, z, e) - R 3 (A, z, e)) ^= exp (0(z" 1+4t ) + 0(A 4e z 2 ^)) (1 + o(l)). 

Written in this form, #3 can be easily compared with i?2, i-e., the formula (|2.4|) . Further we will see 
later that #3 will be dominant if z < (4 — <5) log 1/A for an arbitrarily small constant 5 > 0. Then we 
only have to consider the convergence of H 3 for z < 0(log 1/A). For z < 0(log 1/A), i?i(A, z, e) [ 0, 
R 2 (A, z, e) I and i?3(A, z, e) j as A — > 0. Finally we obtain 

(2.6) ff 3 = ^Lexp(-J)(l + o(l)) if z<0(logl/A). 



Calculation of #4: 

From Lemma 12. 41 

■1-A/e' 



H 4 < A(ip A )exp(-zA(l - A/e'))dip 

J ip A 

< A^-'A^a) exp{-zA(l - A/e')). 
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Here A/A < 1/e'. Then by Lemma f2~Tl and ITU as A -> 



in A/^ ^{e'(i + /?) + !} i n 



where u>(-) is a nonnegative function such as 



lim wit) = 0. 



Finally we obtain 



(2.7) H 4 < 0(A 1 / 2 " e exp(-^(l-u;(e / )) 

For the comparison of H 4 with #3 and Hq, which is needed in later argument, we state some properties 
of H4. For a given small number 5 > 0, there exists e' such that 5 > w(e'). Then for z < 41og 1/A, 

A 5 = exp ^-^ x 41og 1/A^j < exp (^~^ z 

Consequently 

(2.8) Hi < O (A 1 ' 2 - 5 - 6 exp (-|(1 + 5 - w(e')) 
In this representation f/4 can be easily compared with #3. For z > 4 log 1/A 

(2.9) H 4 < 0(Aexp(-cz)). 

In this representation H4 can be easily compared with Hq. 
Calculation of H5: 

From Lemma 12. II for e' < A/ A < 1/e', A{ip) is bounded. Then 

/■l-Ae' 

(2.10) H 5 = A(^)exp(-zA(cp))d(p 

Jl-A/e' 

< (A/e'-Ae) sup A(y>) exp(-M(^)) 

^6(l-A/e',l-Ae') 

< 0(Aexp(-cz)). 

Calculation of H§: 

From Lemma 12. II for A/ A < e' 

M/A\l/(1-A) 

A(l-X;a,P) = { (l + 0( e ') + 0(Alog(l/A))) 

= w 7 +/3 (i+^(A,o), 

where Ri(A,e') = O(e') + 0(A log(l/A)). Substituting the above results into Hq and replacing 
(p = 1 — A by A, we obtain 



8 



Furthermore, replacing 

A = AI 2±<l£ w 

l + R A z 



by x and denning 

/ \ /1/(1-A) A + Af, 
#6 is written as 

A/i i o\i— A rg{z) 

H 6 = { V? / xexp(-x)x^ A dx x {(1 + i? 4 ) (A ~ 1) (l - A)}. 

z Jo 

Then the integration of the right side of the equation is evaluated as 

[ 9{Z) xeM-x)dx + 0(A) = [-xe~ x - e- x ] 9{z) + 0{A) 
Jo 

= l + 0(ge- 9 ) +0(A). 

Here for arbitrarily small constant e 1 and any A — > 0, ^(z) — ► oo as z — > oo. Therefore 

(2.11) H 6 = A(l + /?)z A " 2 (l + o(l)). 

For fixed e and e' the formulas (J23J), ((231), ((213), l(2~7j) . (|2~TT7|l and (f2~TTTl are proved. 

Comparing l(2~3|) . (f^lj) . (|2~5|) . (|2~7|) . (|2~TU)) and (j2~TT|) . the order of H h i = 1, . . . , 6, is summarized 
as follows: 

max (Hi,H 2 ,H4,H 5 ) = o(ma,x(H 3 , H 6 )). 



The relative dominance of H% (|2.6|) and Hq (|2.11|) depends on the value of z. For an arbitrarily small 

5 > 0, 

tf 6 = o(# 3 ) if z < (4 - 5) log 1/A, 
H 3 = o(H 6 ) if z > (4 + 5) log 1/A. 

These are easily confirmed if we substitute (4 ± 5) log 1/A into z. 

Now for proving the uniformity in A and z we let e, e' [ depending on A and z. We first let 
e | 0. We can take e = (log 1/A) 5 " 1 for a constant 5 E (0, 1). Concerning Lemmas 12.112.41 we must 
consider the relations of A 2e > A/A 2 in Lemma 12.31 In the other lemmas the results do not change. 
Substituting e = (log 1/A)" 5-1 into A 2e and taking logarithm of A 2e 

logA 2e = -2(logl/A)' 5 -> -oo. 

Hence A = o(A 2 ) and the conclusion of Lemma 12.31 is not changed. Concerning the calculations of 
Hi, the order of dominance is not affected by the limiting operation e j 0. For the comparison of H$ 
and Hq we may substitute e = (log 1/A) 5 " 1 and z = (4 ± 6) log 1/A into (E^J) and (t2~TTl) . The other 
comparisons are easy. Then we have only to confirm that #3 satisfies (|2.6|) because the calculation 
of Hq is not affected by e. If we substitute e = (log 1/A) 5-1 into (|2,5|) . we obtain ()2.6|) since 

zA 1 - 26 -> 0, /\}/ 2 - t z 1 / 2 + T 0, A 4e z 2r -> 0, for 2 < 4 log 1/A 
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imply 



Ri(A,z,e) | 0, R 2 (A,z,e) j 0, R 3 (A,z,e) | 0, 0(A 4e z 2r ) -> 0, for z < 41og 1/A. 

Next we consider e' | depending on A. We may consider e' = (log 1/A) 5 '" 1 for an arbitrarily 
small 5' G (0, 1). Concerning Lemma 12.11 e ; must satisfy A/e' < A 1 / 2 " 6 . This is easily confirmed, e' 
does not affect the other lemmas. Concerning the calculations H^s, we must consider H^, H§ and 
Hq. The other H^s are not affected by e'. From 1)2.7(1 and the definition of w(e'), the calculation of 
Hi concerning the order, i.e., formulas Ij2.8() and 1)2.9)1 are not changed. We investigate H§ in detail. 
If [1 - A/e, 1 - Ae'), then e' < A/A < 1/e' and A(<p) of Lemma UHl satisfies 

l+fSl^W-*) < AW). 

Hence 

F 5 < (A/e'-Ae')x sup A(</>) exp(-zA(y?)) 

¥?6(l-A/e,l-Ae) 

< O f(A/e' - Ae') x e' 1 ^ 3 exp f-cze' 1 ^ 



< 0(Aexp(-cze' 1/(1 " A) 
Then for z < (4 — 45") log 1/A and for an arbitrarily small positive constant 5" , 

#5 < O(A) < O (A 5 " exp (-£)) = o(i? 3 ). 

For 2 > (4 -45") log 1/A 

# 5 = 0(Aexp(-cz' 5 ')) = o(£T 6 ). 

Hence the order and dominant terms of H^s are not changed. For the calculation of Hq, g{z) must 
go to oo as z — > oo. Considering z > (4 + 5) log 1/A, we can easily confirm this condition. Thus we 
can let e' J, depending on A. Therefore convergence of Hi, i = 1, . . . , 6, only depends on A and z 
and is faster as A — > and z — ► oo. 

Finally, substituting z = (x — Q a /( a— into #3 and Hq and multiplying these by a/(2(a — 
l))(x — C) j we obtain the result for x — £ > 0. For x — £ < we utilize the relations 

/(x; a, /3) = /(— x; a, — /?) of general stable distributions. □ 



Proof of Theorem 2.2 From 1)1.4(1 the derivative of the density is calculated as follows. 

(2.12) f(x;a,(3) = -L-(x - C)" 1 f(x;a, (3) - (x - C) 2/(q - 1} x H', 

a — 1 2(a — l) z 

where 

F' = y ^ 2 (^;a,/3)exp(-(x-Cr /(a - 1) ^;a,/J))^. 

Since behavior of f(x;a,f3) was obtained in Theorem 12. H we evaluate the integration of the second 
term of 1)2.12(1 . As in Theorem 12.11 we divide integral H' into subintegrals, 

k=l 
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where each H' k corresponds to the k-th interval of [0,1 — fi), [1 — fJ>,fo)i [<P0iPa), [<PA,1 — A/e'), 
[1 - A/e', 1 - Ae') and [1 - Ae', 1]. From Lemmas ITT1 l2~2l and IP for A < A/e' 

(2.3) - A) . Wg|±^ (1 + 0(A lo g (l/A))), 



and for v? < v?a ^> A 1 / 2 ^ < A 



^ 2 (^a) = ^ + ^A 1 - 2e + (A), 



A 2 ( l-A) = l + ^A 2 + o( A 2 ). 



Considering A 2 (ip, A) term of the integrand, we easily find the following results on H[ similar to 
H{, i = 1, . . . , 5, in the proof of Theorem 12, II 

(2.14) H[ < 0(exp(- 7 z)/z 2 ), for 7 >i 

(2.15) H' 2 < o(exp(-^ + 0(zA 1 - 2e ) + 0(z 1 / 2+T A 1 / 2 -*)-cz 2T 



(2.16) H> = 1 exp + o(l)), 



(2.17) flj < o(A 1 / 2 -exp(-^l- U ;( e ')))) 

(2.18) < 0(Aexp(-cz)). 



Later we see H^ and H'g dominate the others terms as if 3 and H$ dominate the other H^s in Theorem 
12.11 Since the calculation H' Q is somewhat different, we treat H' & separately. From (|2.13|) . for A/ A < e' 



2 



A 2 (l-X;a,f3) = (1 + ^ 4 (A,0) 2 , 

where f?4(A,e') is defined in the proof of Theorem 12.11 Substituting the above result into H' e and 
replacing ip = 1 — A by A, we obtain 

I±M 2 f Ae ' / WA x2/(l-A) ovn (_J \ + M , WA U/(1-A) 



I (A/A)— exp - z ^, (VA) W-., )(iA . 



Furthermore, replacing 



by rr and defining 



//q is written as 



A = A 



1 + /3 x^" A 



1 + i? 4 ^ 
/1/(1-A) / 1 + -R4 

1 + 



flj = A(1 + /? / A ^ 2 exp(-x)x- A ( i(r x {(1 + ^ 4 ) (A ' 1) (1 - A)}. 

z Jo 
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Then the integration of the right side of the equation is evaluated as 

[ 9 '* x 2 exp(-x)dx + O(A) = [-x 2 e- x - 2xe~ x - 2e- x ] g(z) + 0(A) 
Jo 

= 2 + 0(g 2 e~3) + 0(A). 

Here for arbitrarily small constant e' and any A —> 0, g(z) — > oo as z — ► oo. Therefore 

(2.19) ^ = 2(l + /3)Az A " 3 (l + o(l)). 

Hence for fixed e and e' formulas (|2.14l) - (|2.19|) are proved. 
From fllM - TFT^i 

max(F|,^,^,^) = o(max(^,^)). 
The relative dominance H% and -ff 6 depend on the value of z. For an arbitrarily small 5 > 

tf 6 = o(fl£) if z < (4 - 5) log 1/A, 
H^ = o(H' 6 ) if 2 > (4 + 5) log 1/A. 

These are easily confirmed if we substitute (4 ± 5) log 1/A into z. 

Substituting z = (x — C)"/^ -1 ) into H' 3 and and multiplying each of them by a 2 /(2(a — 
l) 2 )(x — <^) 2 /( a ~ 1 ) we obtain the second integral of the equation l|2.12|) . Combining this with Theorem 
12.11 we obtain the desired result. For x — Q < we use the relation f(x;a,(3) = f(—x;a,—j3). For 
proving uniformity in A and z we need to let e, e' [ 0. We omit the derivation because it is quite 
similar to that in the proof of Theorem 12.11 
□ 

From Theorem l2.ll Theorem l2.2l and formulas 1)1. 2J) . 1)1-3)1 . we obtain behavior of the score functions as 
a | 2. These score functions in the corollaries below are needed for obtaining the Fisher information 
matrix of the general stable distributions. 

Corollary 2.1 Under the same conditions and notations of Theorem \2.1\ and \2. 21 define 

g^{x;a,f3) = -g'(x;a,(3)/g(x;a,(3). 
Then for an arbitrarily small e > 0, there exist Aq and xq such that for all A < Aq and \x\ > xo, 

\fJx;a,(3)/gJx;a,(3) - 1| < e. 



Furthermore define 



F 2 (x;a,f3) 



2 : 
3 



Then 

g^(x;a,P) 
g(x;a,f3) 



Ffa;a,{3)(l + o(A 5 / 2 )) if \x - C| < (2 - <5)(log 1/A) 1 / 2 
F 2 (x; a, /5)(1 + o(A 5 / 2 )) if \x - (\ > (2 + 5) (log 1/A) 1 / 2 , 
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uniformly in \x\ > xq. For the remaining interval (2 - J) (log 1/A) 1 / 2 < x < (2 + 5) (log 1/A) 1 / 2 , 



< cx 



g(x;a,f3) 
uniformly in \x\ > a^o- 

Corollary 2.2 Under the same conditions and notations of Theorems \2.1\ and \2.2i define 

g a (x; a, f3) = -g(x; a, 0) - (x - ()g'(x; a, f3). 
Then for an arbitrarily small d > 0, there exist Aq and xq such that for all A < Ao and \x\ > xq 

\f a (x;a,f3)/g a (x;a,f3) - 1| < e. 



Furthermore define 



F'(x;a,(3) 
F*(x;a,j3) 



(x - ()< 



Then 



g a (x;a,f3) 
g(x;a,(3) 



F (T 1 (x;a,/3)(l+o(A 1 /2)) if | x - C| < (2 - 5)(log 1/A) 1 / 2 
F CT 2 (x;a,/3)(l+o(A 1 /2)) if | x - C| > (2 + 5)(log 1/A) 1 / 2 

For the remaining interval (2 — 5) (log 1/A) 1 / 2 < x < (2 + 5) (log 1/A) 1 / 2 , 

g a (x;a,/3) 



g(x;a,P) 



<cx 2 



uniformly in \x\ > xq. 



3 Derivatives of density with respect to the parameters a and j3 

In this section we obtain the derivatives of the density with respect to the parameters a and j3 by 
analyzing the inversion formula. The density of general stable distributions can be written as 



(3.1) 



i r°° 

f(x; a, (3) = -Re / e- ltx $(t; a, (3)dt, 
n Jo 



where $(i;a,/3) is defined in f|l . 1|) . Utilizing ()3.1|) . we derive the following two lemmas concerning 
derivatives with respect to a and (3, which are extensions of Lemma 5 of Nagaev and Shkol'nik (1988) 
to non-symmetric case. 

Lemma 3.1 As A = 2 — a — > 0, there exists xq and for all \x\ > xq, 

M 3 M 4 + M 5 log \y\ 



f a (x;a,P) 



\y 



^{l + / 9* + A(M 1 + M 2 log|y|) + ^ + 



\y\ 



where y = x — £ and M\ , M$ are some constants. 
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Lemma 3.2 Under the same conditions and notations of Lemma UJ.ll 

Asgny / M 7 M 8 

^(x; a, /?) = ^1 + AM 6 + — + — 

where Mq, M7 and Ms are some constants. Furthermore 

\ffs(x;a,P)\=0(A), 

uniformly x £ R. 

Proof of Lemma 3.1 First we consider the case of y = x — £ > 0. Differentiating / by a, we get 
(3.2) f a (x;a,(3) = ^ J°° e~ ity exp (-t a (l - z/3tan (^))) 

x {-nogt (l - */3tan )) + (f - t)^/ (cos (^)) 2 } dt. 
Transforming f a (x;a,P) by t — ► s/y, we obtain 
f«(x;a,P) = ^Re^ 0O e - is exp(- S a y- Q (l-i/3tan(^)) 2 ) 

x |- a «y-«(loga + logy) (l - i/3tan (^)) + ( S V° - sy" 1 )^/ (cos (^)) 2 } ds 

Further transforming f a (x; a, (5) by s — > te l(; , c = — ^- as in Lemma 5 of Nagaev and Shkol'nik (1988), 
we find 

f a (x;a,f3) = ^Ree i? | exp (-^e i? + it a y~ a (l - iptun (^))) 

x |^y- a (logt + k + logy) (l - t/Stan (^)) - (it a y- a + ty^e 1 ^/ (cos ^)'| 
= K x + K 2t 

where K\ and are obtained in formulas (|E,2|) and (jE.4|) of appendixEl Then the lemma is proved 

for y = x — C > 0. For y = x — £ < 0, we use the relation / a (x; a, /?) = f a (—x; a, — /?). 

□ 



Proof of Lemma 3.2 First we assume y = x — ( > 0. Differentiating the density Ij.S.lj) by (5 and 
doing the same transform as in Lemma 13.11 we obtain 

a, 0) = f 4 Re r e*e-**t*dt - -Re f° if*e**tdt + + . 

vry Vy jo y Jo y y 1+c v 

The calculation of fp(x;a,P) is easy by the equations (|E.1|) and (|E.3|) in appendix IE1 As A — > 0, 
tan ^y- = — f A + o(A 2 ). Finally we utilize the relation fp(x; a, (3) = —fp(x; a, — /?) and the first part 
of the lemma is proved. The second part of the lemma is obvious from 

I tan — I f'°° 

\fp(x;a,P)\< 1 - 2i / e~ ta \t a -t\dt. 

n Jo 

□ 
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Table 1: Limit of information matrix at a = 2 



he 




a 


a 


P 




0.5 











a 





2.0 


— oo 





a 





— oo 


DO 





P 















4 Information matrix of general stable distributions 



The Fisher information matrix of general stable distributions are derived in this section. Table^gives 
limiting values of the information matrix at a = 2. To the author's knowledge other than the diagonal 
element and 1^, these limiting results have not been given in literature. In following theorem we 
obtain asymptotic behavior of information matrix in more detail for the important diverging cases. 

Theorem 4.1 As A = 2 — a — > 0, behavior of the Fisher information matrix of general stable 
distributions at \i = 0, a = 1, (3 ^ ±1, is given as follows. 







h a 


hp 


* 


ha 


ha 


h/3 


* 


* 


ha 


hf3 


* 


* 


* 


hp . 



0.5 + o(l) o(l) 
* 2.0 + o(l) 



oil) 



O(A) 



- 1 log log 1/ A o(Aloglogl/A) 
x(l + o(l)) 



4A logl/A 

x(l + o(l)) 



log 1/A 



4(1-^) logl/A 

Note that for and the rates of convergences to are not obtained and only the limiting 
values are obtained. For I aa , I aa and hp the exact rates of convergences are obtained. However for 



symmetric stable distributions since i. 



\J(7 



/.la 



0, the above 3x3 matrix presents exact limiting 
behavior of the Fisher information matrix as A — > 0. The idea of the proof is based on the proof of 
Theorem 2 of Nagaev and ShkoPnik (1988). 

Proof of Theorem 4.1 In our proof we use the following notations. T is a sufficiently large constant 
such that Theorem 12 . II and Corollaries 12.11 and 12.21 are applicable. We also denote 



Xl (A) = (2 - 5) (log 1/A) 1 / 2 , x 2 (A) = (2 + J) (log 1/A) 1 / 2 , x 3 (A) = exp(A 



-1/2n 



where 5 > is an arbitrarily small constant. Further we use the notation "const" for some proper 
constants. We prove Iaa in some detail and the proof of the other Iqq are given in Appendix |FJ . 

Proof of I aa : 

In our proof we utilize Theorem 12. II and Lemma 13. II We divide integral of Fisher information matrix 
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into two parts, 



^ = r {fa(x + C;a,(3)} 2 dx | r™ {f a (x-C;a,-P)} 2 dx 
Jo f(x + (;ot,p) Jo f(x- C,a,-(3) 

= I 1 +I 2 

This is obtained by the relation f(x;a,(3) = f(—x;a,—(3) for x — £ < 0. l\ a is calculated first. 
Further we divide integration I aa into five subintegrals, 

5 

k=l 

where each I aa (k) corresponds to the fc-th interval of [0, T), [T,xi(A)), [xi(A), 22(A)), [22(A), 2:3 (A)) 
and [23(A), oo). 
Clearly 

(4.1) J Qa (l) < oo. 

For I aa {2) 



f(x + (;a,/3) = f(x;2)(l + o(l)), f a (x + (;a,/3) = const x x 

Then 

fa:i(A) /t-2\ i 

OA_fi / X \ . 1 



A-3 



Ai^Clogl/AJB/z-A- 



yxi(Aj f x 2 \ 

(4.2) I aQ (2) = const x y x 2A " 6 exp ( — J cfe < const x 

For I aa (3) 

f(x + (;(3,a) = {(1 + /3)Ax A " 3 + /(x; 2)}(1 + o(l)), / Q (x + C; a, /?) = const x x A " 3 . 

Then 

(4.3) 7 aQ ,(3) < const x — - / x dx = const x 

A J XT (A) 

For J aa (4) 



1 T 2(A) A_ a , . (5 

i Xl (A) Alogl/A" 



/(x + C; a, /?) = (1 + /?)Ax A " 3 (l + o(l)), / a (x + C; a, /?) = -(1 + /?)x A " 3 (l + o(l)). 

Then 



(4.4) 7 aa (4) = — ^ / x^dx(l + o(l)) = — " (1 + o(l)). 

A 7 X2 ( A ) 8Alogl/A 

For 7 aa (5) 

f(x + C; a, /?) = (! + (3)Ax A ~ 3 (l + o(l)), / Q (x + C; a,/J) = -(1 + (3 + Alog x)x A ~ 3 (l + o(l)). 
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Then 



1 f°° 

(4.5) I aa (5) = const x — / £ A " 3 {max(l + (3, Alogx)} 2 dx 

A Jx 3 (A) 

„ e (l+/3)/A 1 

< const x A / x A ~ 3 (log x) 2 dx + const x — / 

Jx-a(A) A Jet 1 



1 Z" 00 

< const x — / x A_3 dx 

= 0(e~ 2 / Al/2 /A) -^0, as A -»■ 0. 
From the formulas (|4.1j) - ([4.5|) we obtain 



x dx 

23(A) ^ Je^+P)/* 

oo 



Setting — > —(3 in l\ a , we obtain I 2 a . Adding l\ a and I 2 a , we prove the assertion. 
□ 

A Proof of Lemma 12.11 

As A — ► 0, ag = ^ arctan (/3 tan(^)) in A(<^, a, (3) is expanded as 

ag = l arctan (/3tan (^)) = -/3A - = ^ (if A 3 + o(A 4 ). 

Since the smaller remainder terms other than — (3A are not needed in the following argument we 
substitute — (3 A into A(tp, a, f3) for ag in advance. Then for convenience we write A(ip, a, (3) of ()1.5|) 
as 

cos(f {(1 - A)y> - gA|) / cos(f <p) cos(f /?A) \ V^ 1 "^ 
1 • j ^' ^ j sin(f {(2 - A ^ - /3A}) V,sm(f {(2 - A)<p - (3 A}) J 

C(<p,A,f3) ( D(ip)E(A,(3) \ 1/{1 - A) 
B(<p,A,/3) V B(<p,A,/3) J 



where 



B(<p,A } 0) = sin(|{(2-AV-/3A} 
C(<p,A,f3) = cos(|{(l-A)^-/3A}) 
D(p) = cos (?-(p) , 



E(A,(3) = cosf^A 
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Then substituting 1 — A for ip we obtain 

(A.2) B(1-\,A,0) = sinQ(A + 2A- AA + /3A) 

(A.3) C(l-\,A,/3) = sin(|(A + A- AA + /3A)) 

(A.4) D(l-\) = sin(|A), 

(A.5) E(A,(3) = cos(|/3A 

Note that if A — > 0, then A — > since < A < A/e' and e' > is an arbitrarily small fixed number 
(or even when e' converges to slower than A). Expanding in terms of A — > 0, we have 

B(1-X,A,(3) = ^a(i + 2^-X + 0\ +0(A 3 ), 

C(1-X,A,P) = |A^1 + ^-* + P) +0(A 3 ), 

D(l-A) = |AxA + 0(A 3), 

£(A,/3) = l-^g/3A) 2 + o(A 3 ). 

Substituting these expansions to A(<p, A, /3) of (|A.lj) and utilizing xVCA-i) = x~ 1 — x~ l logxA + o(A), 

we prove the desired result. 

□ 

B Proof of Lemma 12.21 

We prove A((p A , A,/3) first. When ip = ip A = 1 — A 1 / 2 ^, we replace A by A 1 / 2 " 6 in the formulas 
□a- CO)> . Then 

B(<p A ,A,0) = sing{2A 1 / 2 - + (l + /3)A-A 3 / 2 -}), 
C(<p A ,A,0) = sing{A 1 / 2 - + (i + /?)A-A 3 / 2 - e } 
D( VA ) = singA 1 / 2 - 

Note that E(A, (3) is the same as in Proof of Lemma 12. II and we omit it. As A — ► we have 
(B.l) B(p A ,A,p) 

= vrA 1 / 2 - |l + 1(1 + /3)Al +e - y A1 ^ -\ A ~ T (1 + ^ A3/2 " € + °( a3/2 )} . 
(B.2) C(<p A ,A,p) 

= \A^ |l + (1 + ^)A 1 / 2+e - l^A 1 ^ - A - y (1 + /3)A 3 / 2 ^ + o(A 3 / 2 )} , 

(B.3) i% A ) = |AV 2 - (l - I^A 1 ^ + o(A 3 / 2 )) . 
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Substituting these expansions into A(<p,A,(3) of (TOll . we obtain the desired result. 

Next we prove assertion concerning A'((p A , A, (3). From the equation ([1.5)1 . we can write 

(B.4)A / ( v ,A ) /3)=A(^,A,/3) 

7T f (2- A) 2 cos(f{(2- A)vj-/3A}) sin(|{(l - A)cp — ft A}) 1 sin(f^) 



2{ A-l sin(f {(2 - A)p - /3A}) ^ A ^cos(f {(1 - A)ip - pA}) l-Acos(f^) 
^' A)/3) X 2 1 A-l B((p,A,/3) ~ (1 " A) CVA^) " T3A%) j ' 



where 



F(<p,A,0) = cos(|{(2- A)v9-/3A}) , 
G(^,A,/3) = sin(|{(l-AV-/3A}), 



£%) = sin 



Then substituting tp A for c/?, we obtain 



F(<p A ,A,p) = -cos(-{2A 1 / 2 - + (l + /3)A-A 3 / 2 -} 
G(<p*,A,p) = cosglA^ + OL + ^A-A 3 / 2 -} 
H(^a) = cosf^A 1 / 2 ' 



As A -► 0, 

(B.5) F(<pa, A, P) = -l + I 7 r 2 A 1 - 2e + i 7 r 2 (l + /3)A 3 / 2 - e + (A 3 / 2 ), 

(B.6) G((fA, A, P) = l-^A 1 - 2e -^(l + /3)A 3 / 2 - £ + o(A 3 / 2 ), 

2 

(B.7) f% A ) = l-^A 1 ~ 2e + (A 3 / 2 ). 



Combining l(B~l |) -(|B~3 j) , (|B~5|) - (|B~T|) . A(<p, A,/3) and A'(ip,a,P) of (|B~4)) . we prove the lemma. 
□ 

C Proof of Lemma 12.31 

We further expand B(<p,A,P), C((p,A,P) and D(<p,A) in (|A~2 )) -(|A~4 jl with respect to A. As A -> 
and A ^ such that A 1 / 2-6 < A, 



(C.l) B(l-A,A,/3) 
(C.2) C(l-A,A,/3) 



2 I 6 V2 / 120 \2 / A 
(C.3) £>(1-A,A) 
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Substituting these equations to A((p, A, /3), we obtain the result. 
Concerning A'((p,a,j3) we use formulas 

F(1-X,A,P) = -cos(^{2A + (l- A + /3)A} 
G(l-A,A,/3) = cos (|{A + (1 - A + p)A} 
H(l-X) = cosgA). 

As A -> and A -> under A 1 / 2 " 6 < A, 

(C.4) F(l-A,A,/3) = -l + I(7rA) 2 + o(A 3 ), 

(C.5) G(l-A,A,/3) = i_IgA) 2 + o(A 3 ), 

(C.6) F(l-A) = l-IgA) 2 + o(A 3 ). 

Substituting (|C.1 |) -()C.3 | ) and l|C.4|) -()C.6 | ) to A'(<p,a,f3) of l)B.4|) . we obtain the desired assertion. 
From the equation (|1.5jl or A'(p, a, 0) of (|B.4|) . we can write 

(C.7)A"(^A,/3) 

-%^)x!( ( 2 ; A fS-^! -(i-A)^^) 1 



2[ A-l B(<p,A,j3) v 'C(<p,A,P) 1-A%) 
vr 2 f (2 -A) 3 (1-A) 3 



4(1 - A) [ ( S in(f {(2 - A)<^ - /3A})) 2 (cos(f {(1 - A)y> - /?A})) 2 (cos(f 

a r vr f (2-A) 2 F(y,A,/3) A;> G(y,A,/?) 1 fffc) 1 

A ^ A ^ X 2\ A-l ^,A,/3)- (1 - A) ^A^)-r3A^)j 
vr 2 f (2 -A) 3 (1-A) 3 1 



+A&) 



4(1 - A) 1 £? 2 (^, A, /?) C%, A, /3) -D 2 (y) 



Combining <f(TT |t - (fTTa jl . ft ^ -ftT fy . A'(tp,a,P) of ((B3J) and A"(ip,a,p) of (EH), we obtain the 

desired assertion. 

□ 



D Proof of Lemma 12.41 

We only need to prove for x — C > 0, because for x — £ < /(x; a, /3) = f(—x; a, —p) and £* in 1)1. 5j) 
does not change, namely 

to / /7ra\\ . , i it a ( /ira 

g = sgn(a; — Q)-^— arctan I ptan I — I I = sgn(— x + Cj-^r- arctan I — p tan I — 



We use notation Pb G (—1, 1) which corresponds to the parameter P E (—1, 1) for (B) representation 
of characteristic functions. The reason is that the proof in terms of (B) representation is simpler since 
ag = ~PbA (see p. 74 of Zolotarev 1986 ; the notation g in our paper corresponds to 9 in Zolotarev) 
and we do not need to expand ag with respect to P as in the proof of Lemma 12.31 Moreover our 
purpose is only to determine the sign of A'(tp). We divide the integral interval as (— g, 1 — A/e'], 
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[1 — A/e', 1) and prove the lemma for each interval. Except for the interval (— g, 0) we follow the line 
of the proof of Lemma 4 of Nagaev and Shkol'nik (1988). For (1 — A/e', 1) the proof is obvious from 
Lemma 12. II For (— g, 1 — A/e'] we have to consider the sign of g. Note that 

sgn(p) = -sgn(/3 B ). 

First for < -g < if < 1 - A/e' ^ A/e' < A < 1 + g, we get 

where 

Ai(l-A) = -(2- A) 2 cot (|(A + 2A- AA + /? B A) 

+ (1 - A) 2 cot (|(A + A - AA + £ B A)) + cot (|a 

From the formula (4.3.70) in [4] we have 

Ah x , 2(l + fe) 2 A 2 ~ 

" A) = 7rA(A + 2A - AA + (3 B A)(A +X-AX+ (3 B A) + ^ 47fc(A) ' 

where c^'s are positive numbers and 

7fc(A) = (^a) 2 " +1 [(2 - A) 2 {2 - A + (1 + fi B )A/X} - (1 - A) 2 {1 - A + (1 + f3 B )A/X}] . 

Since A/A < e' for arbitrarily small e' > 0, A\{ip) > is easily confirmed. 

Secondly we investigate the case — g < ip < 44> ay? = 7A where 7 £ (/3s, 0) for f3 B < 0. Then 

Ai(tp) = (2- A) 2 cot (|a((^ + £.)) +(l-A) 2 tan(|{a£.+ (a-l)v9}) +tan(| e 

= (2 - A) 2 cot fj( 7 - (3 B )A) + (1 - A) 2 tan f J ( 7 - ft, - + tan ^ 



.2 W ' ' / v y V 2 V ' 2- A/ / V 2 2 

As A 0, 

(2 - A) 2 

We easily find that A\{(p) > for sufficiently small A. 
□ 

E Calculations of f a 

For all y > y$ = xq — (, we have 



exp ( ^ a y" Q ( 1 - i/3 tan (^— J J J = 1 + c t a y~ a (1 - i(3 tan ^ 
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where cq G C is some complex constant. Utilizing this expansion in f a (x;a,(3) we obtain 



Jfi = -Re e t? exp -iie K + l-i/?tan. 

vry 7 V V V 2 / 

xit Q y"°(logt + k + logy) (l - i/3tan (— )) ^ 

= -L^Re e- e"^ {l + c fy- (l - z/3tan (^)) } 

xit a (logt + K + logy) (l -i/3tan (^)) dt 

= -4— Re e i? /"°° e" ite " it Q (log t + + log y) f 1 - t/3 tan ( ) dt 
ny 1+a Jo V \ 2 J J 

l^Re e* I" e-^ico t 2a (l - /5tan (^))" (logt + k + log y)dt 

J 



Try 1 

1 „ , v /"°° s, A(c + clogy) c + clogy 

—— Re e K / e~ lte it a log i + + log z)dt + — * -— ^± + rr^- 

vry 1+Q Jo y 1+a y 1+2a 

Note that As A -> 0, tan 3p = _| A + o(A 2 ). Then we calculate the integrations in i^i above. Write 
Ai = cos? and A2 = sin?. Here at a = 2, Ai = l/y/2 and A2 = —l/y/2. Prom the formulas (3.5), 
(4.40) in part I and (3.7), (4.17) in part II of Oberhettinger (1990), we have 

poc poo 

Re e !? / e- ite "it a log tdt = / e tx 'H 2 log i{Ai sin(iAi) - A 2 cos(tAi)}dt + Ac 
J ■/ 

= 4AiA 2 arctan + 8AiA 2 (A| - A?) Q - EU^) + Ac 

H A 

= -5 + Ac, 

f'OO f'OO 

(E.l) -Re?e i? / e~ itel 't a dt = - e tx H 2 {\ 2 sin(iAi) + Ai cos(tAi)}d* + Ac 

= -4?AiA 2 + Ac 
= -2+ Ac, 

/*oo /*oo 

ReeM e- ite% \t a dt = / e iA2 t 2 {Ai sin(tAi) - A 2 cos(tAi)}dt + Ac 

= -2Af + 2A| + Ac 
= Ac, 

where EU means Euler's constant. Substituting these equations into K±, we obtain 
(E.2) Kl = - JL { 1 + A(c + c log y) + C -±^ 
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For K2 we do the similar calculations. 



x(A ° + ty- 1 e«)^-/(oofl T ) ^ 

P_ t-> _ _ic /" _— iie* c f 1 1 — ci /1 .'oi._ fKOt 

2y(cos(^))2 



rRe e if e~ ite " {l + ct Q y" Q (l - i/3tan (™)) } (* a 2T a - ity^e^dt 



I —Re /°° e~ ite *t a dt - -Re e 2 ^ /°° e~ ite *itdt + + ° 



,i+ft 



2y(cos(^)) 2 U a y 2/ 7o y 2a 2/ 

Substituting the following two equations into i^j we obtain i^- The first equation is obtained from 
formula (|E.1|) in the calculation of i^i. The second equation is obtained from the formulas (3.5) in 
part I and (3.7) in part II of Oberhettinger (1990). 

/•oo 

Re e k / e~ iteK t a dt = -2 + Ac, 
Jo 

POO POO 

(E.3) -Reie 2i " e~ iie " ' tdt = / e* A2 t{2AiA 2 cos(tAi) + (A 2 , - Xj) sin{t\i)}dt 

Jo Jo 

= 2A 1 A 2 (A 2 -A 2 ) + (A 2 -A 2 )(-2A 1 A 2 ) 
= 0. 

Hence 

(E.4) ^ 2 = ^(_l + Ac +4 + - 

yl+a I yft y 

F Proof of I 00 other than I aa 

Here we give proofs for IppJapJaaJap, Ifi/3 and I m . In the following we assume T > xq where xq is 
defined in Theorem 12.11 and Corollaries 12.11 and 12.21 

Proof of Ipp : 

Here we utilize Theorem 12 . 1 1 and Lemma 13.21 However the proof of Ipp is essentially the same as the 
proof of I aa because 

fp(x;a,(3) = -j^-sga(x - Qf a (x; a, (3){l + o(l)) 

for sufficiently large \x — (\. Therefore we omit the detailed derivation of Ipp. 
Proof of I a p : 

Here we utilize Theorem 12.11 Lemma 13. II and Lemma 13.21 From the same reason as Ipp we omit the 
proof. 

Proof of I aa : 

We make use of Lemma 13. II and Corollary 12.21 We divide integration I aa into two subintegrals, 
f a {x + (;a,(3)f a (x + (;a,P) dx+ f°° f a (x - C; a, -P)f a (x - Cj a, -@) dx 



f(x + (;a,(3) J f(x - (;a,-f3) 

I 1 +I 2 

^aa 1 (jot 
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is calculated first. Further we divide integration l\ into four subintegrals, 



I act ~ ^ j lgajk), 



k=l 



where each I aa {k) corresponds to the integration of l\ a for the fe-th interval of [0, T), [T, x\(A)), 
[xi(A), X2(A)) and [a;2(A),oo). For 1^(1) clearly 



(F.l) 

For I aa (2) 

Then 
(F.2) 

For I aa (3) 
and 



1^(1) < oo. 



-(1+/?) 
1 + /? 



x(A) X A-1 



dx(l + o(l)) 



rp 2 

loglogl/A(l + o(l)). 



f a (x + C,«,/3) 



< const x x 2 , \f a (x + (;a,(3)\ = (1 + (3)x A ~ 3 {l + o(l)), 

const x x A_1 (l + o(l)). 



/(x + C,«,/3) 

+ C; a, /3)/a0 + C; a, /?) 



/(x + C;a,/?) 



It is easy to see 
(F.3) 

For I aa {4) 



I aa (3) = const x / x A ~ l dx = 0{5). 
Jxi(A) 



f a (x + C,a,(3) 



<2(l + o(l)), \f a (x + (;a,P)\ = (1 + /3)x A " 3 (l + o(l)). 



f(x + (;a,P) 
Thus we easily show 

(F.4) W4)-0, asA^O. 

From (JEU), (lF^|) . (lF~3l) and (ITU) , we obtain 

£ a = ~(l + /3)loglogl/A(l + o(l)). 

Setting /? — > —(3 in we obtain I 2 a . Adding l\ a and I 2 a , we prove the assertion. 
Proof of I a p : 
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Here we utilize Lemma 13.21 and Corollary 12.21 From the same reason as Ipn we omit the proof. Note 
that the result can be expected from behavior of A x I aa and the asymmetry of fp{x;a,P) around 
C, which we see in Lemma 13.21 



Proof of I. 



Here we utilize Lemma 13.21 and Corollary 12.11 We divide integration I„n into two subintegrals. 



'«3 



f^x + C; a, /3)fp(x + C; a, P) 



dx + 



ffjjx - C; a, -/3)fp(x - C; a, ~P) 
fix - (;a,-P) 



dx 



T 1 4- T 2 
I Hf3 + 1 w3- 



I^p is calculated first. Further we divide integration into two subintegrals, 



w3 



fc=i 



where each Ifj,p(k) corresponds to the integration of 1^ for the k-th interval of [0, T) and [T, oo). For 

f p (x + C;a,p) = 0(A), /^ + C;«,/3) < const 

Then 1^(1) = O(A). For 1^(2) it follows that 

U(x + C;a,/3) 



f{x + (;a,P) 



< const x x, fp(x + C; a, /?) = const x A: 



„A-3 



uniformly in x G [T, oo). Thus we find 1^(2) = 0(A) easily. Setting (3 — > —(3 in 1^, we obtain i^g. 
Adding P~ Q and 7^, we obtain the desired result. 



V/3 

Proof of 7, 



Here we utilize Lemma l3.1l and Corollary 12.11 We divide integration I aa into two subintegrals, 



xl. 



fia 



U(x + (;a,P)f a {x + (;a,P) 



f(x + (;a,P) 



dx 



ffjjx - C; a, -P)f a (x - C; a, -/?) 
fix - (;a,-P) 



dx 



I 1 +/ 2 

fia 1 fia 



lh a is calculated first. From Corollary 12.11 and Lemma 13.21 we can find an integrable function gi(x) 
such that for x € (0, oo) and all sufficiently small A 

U(x + (;a,P)f a (x + (;a,P) 



< 9\{x). 

f(x;a,p) 

Hence we can apply the dominated convergence theorem to I* . Then 

/ M (x;2)/ a (x;2) 



lim I}, 



A^O 



fia 



f(x;2) 



-dx. 



We apply the same arguments as /L, to I 2 and obtain 



lim I 2 



A->0 



fia 



f^x;2)f a (x;2) 
f(x;2) 



dx. 



Adding /L, and I 2 a , we prove our assertion. 



□ 
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